Phonemic pairing glue in cuprates and related high-temperature superconductors 
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Along with some other researches we have realised that the true origin of high-temperature super- 
conductivity should be found in the strong Coulomb repulsion combined with a significant electron- 
phonon interaction. Both interactions are strong (on the order of 1 eV) compared with the low 
Fermi energy of doped carries which makes the conventional BCS-Eliashberg theory inapplicable 
in cuprates and related doped insulators. Based on our recent analytical and numerical results I 
argue that high-temperature superconductivity from repulsion is impossible for any strength of the 
Coulomb interaction. Major steps of our alternative polaronic theory are outlined starting from 
the generic Hamiltonian with the unscreened (bare) Coulomb and electron-phonon interactions 
accounting for critical temperatures of high-temperature superconductors without any adjustable 
parameters. 
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I. HIGH-TEMPERATURE 
SUPERCONDUCTIVITY FROM REPULSION: 
IMPOSSIBILITY 

With a few exceptions [l[ many researchers share the 
view that the electron-phonon interaction (EPI) is insuf- 
ficient for high temperature superconductivity at least 
in the framework of the conventional Bardeen-Cooper- 
Schrieffer (BCS) theory or its intermediate-coupling 
Eliashberg extension ||. Phenomenologically, the pair- 
ing mechanism of carriers could be not only phonemic 
as in the BCS theory or its strong-coupling bipolaronic 
extension [3], but also excitonic, plasmonic, magnetic, ki- 
netic, or due to some purely repulsive interactions com- 
bined with the unconventional pairing symmetry of the 
order parameter. 

Quite naturally, a number of authors assumed that 
the electron-electron interaction in high-temperature su- 
perconductors was strong but repulsive providing high 
T c without phonons via superexchange and/or spin- 
fluctuations in unconventional (e.g. d-wave) pairing 
channel. A motivation for this concept can be found 
in the earlier work by Kohn and Luttinger (KL) 
who showed that the Cooper pairing of fermions with 
any weak repulsion was possible since the two-particle 
interaction induced by many-body effects is attractive 
for pairs with large orbital momenta, 1 » 1. While 
the KL work did not provide the specification of the 
actual angular momentum of condensed Cooper pairs, 
Fay and Layzer Q found that a system of hard-sphere 
fermions condenses at low densities into p-orbital state 
(I = 1). The critical transition temperature T c of re- 
pulsive fermions was estimated well below 0.1 K with 
very little enhancement due to flatness of the Fermi sur- 
face [f|. In two dimensions (2D) the KL effect is absent 



for the parabolic band dispersion, but the d-wave low- 
T c pairing was found with the repulsive Hubbard U (i.e. 
hard-sphere) potential when tight-binding corrections to 
the electron energy spectrum are taken into account [7[. 

More recent studies claimed that even a weak repul- 
sive Hubbard U combined with lattice induced band- 
structure effects results in higher values of T c in the spin 
singlet d-wave channel near half- filling "encouragingly 
similar to what is found in the cuprate high-temperature 
superconductors" [1, Q. 

Different from these studies we have analysed the KL 
problem with the realistic Coulomb repulsion rather than 
with the hard-core Hubbard U [10]. The second-order 
diagrams, responsible for the unconventional pairing in 
the KL theory Fig. (JT|) , are prohibitively difficult to eval- 
uate when the Fourier transform of the repulsive poten- 
tial, v(q) depends on q, so that most previous and recent 
studies [a Q confined to the hard-sphere Hubbard U re- 
pulsion with t>(q)=constant. In that case the first order 
does not contribute to the pairing vertex in any uncon- 
ventional channel with 1^1, and the only (attractive) 
contribution comes from the second-order exchange dia- 
gram providing p-wave pairing in 3D for parabolic dis- 
persion [5| and d-wave pairing in 2D for the tight-binding 
dispersion [7H9(. 

Actually the screening length in the dense Coulomb 
gas, where the perturbation expansion in powers of s = 
e 2 /irhvp makes sense, is large compared with the charac- 
teristic electron wavelength ~ 1/fc.F, so that it is unrea- 
sonable to treat the weak Coulomb repulsion as short- 
ranged as in Ref.Q ( Vp and kp are the Fermi speed 
and the wave- vector, respectively). We have included 
all buble diagrams in the screened Coulomb potential, 
v(q) = Aire 2 /[q 2 + s(q)] with the static Lindhard func- 
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FIG. 1: (Color online) Angular momentum of the order pa- 
rameter as the function of the dimensionless repulsion s = 
e 2 /nhvF in the 3D Coulomb gas. Inset shows the first or- 
der (a) and the second-order (b,c,d,e) contributions to the 
two-particle vertex. (Reproduced from Ref. [10(. ©American 
Physical Society, 2011.) 
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and carefully calculated all second-order diagrams in the 
3D Cooper channel [Ioj |. Then solving the linearised BCS 
equation allowed us to find the symmetry (i.e. the angu- 
lar momentum I) of the order parameter, shown in Fig.([T]) 
as the function of the interaction strength s. Remarkably 
we did not observe p and d-wave pairing in the whole re- 
gion of the Kohn-Luttinger perturbation expansion and 
beyond (up to s = 3). There is a pairing in higher mo- 
mentum states, I ^ 3, but the corresponding eigenvalues 
are numerically so small (A3 « 0.0011 at s = 3), that the 
corresponding T c « Ep exp(— 1/A) is virtually zero for 
any realistic Fermi energy Ep. 

To analyse the unconventional pairing in a two- 
dimensional (2D) electron gas on the square lattice with 
the tight-binding energy dispersion, 

E p = —2t{cos(p x a/h) + cos(p y a/h)} — [i (2) 

(a is the lattice constant) we used the model 2D Coulomb 
potential v(q) = 2ire 2 /(q + k) with a constant inverse 
screening length k. For the half- filled band the Fermi 
level, /i is found at the van-Hove singularity (vHs) of the 
density of states, (J, = 0, so that one might expect a strong 
enhancement of the unconventional T c near half-filling 
due to vHs proximity [H, 0|. Numerically solving the 
linearised BCS equation by discretization of the Fermi 



surface we reproduced fairly well the results of Ref. [1] for 
the 2D Hubbard model, where the ground state is Bi g 
spin singlet with the d-wave symmetry x 2 — y 2 close to 
the half-filling. However, using the screened Coulomb re- 
pulsion instead of the Hubbard one qualitatively changed 
the ground state, so that contrary to Ref. [7H9| neither 
p- nor d-wave pai ring is possible at any filling while the 
repulsion is weak [10|. This surprising result is due to a 
nonvanishing first-order repulsive contribution to uncon- 
ventional channels from the finite range interaction. 

In cuprate and related superconductors the Coulomb 
repulsion is believed rather strong s S> 1, so that the 
perturbative KL approach might have no direct relevance 
to these materials. Different numerical techniques have 
been applied to elucidate the ground state of the repulsive 
Hubbard model in the intermediate to strong-coupling 
regime, sometimes with conflicting conclusions. In par- 
ticular, recent studies by Aimi and Imada [ll| using a 
sign-problem-free Gaussian-basis Monte Carlo (GBMC) 
algorithm showed that the simplest Hubbard model with 
the nearest-neighbor hopping has no superconducting 
condensation energy at optimum doping. This striking 
result has been confirmed in the variational Monte Carlo 
(VMC) studies by Baeriswyl et al. [Hj], who found, how- 
ever, some condensation energy away from optimum dop- 
ing and also adding next-nearest neighbor hopping. 

Using similar VMC simulations of the Hubbard model 
but including a reaslistic finite-range electron-phonon in- 
teraction, we found that even a relatively weak EPI with 
the BCS coupling constant A ps 0.1 induces the d-wave 
superconducting state in the model with the condensa- 
tion energy several times larger than can be obtained 
with the Hubbard repulsion alone [l3[, Fig.©. More- 
over, the unconventional superconductivity has been 
shown to exist due to a finite-range EPI fl3l.[l4| without 
the need for additional mechanisms such as spin fluctua- 
tions. 

Based on these findings one concludes that the p- and 
d-wave Cooper pairing from the weak Coulomb repul- 
sion is not possible between fermions in any dimension. 
Pairing in higher momentum states (I ^ 3) has the corre- 
sponding T c virtually zero for any realistic Fermi energy. 
The unconventional pairing from the strong Coulomb re- 
pulsion is highly unrealistic either since the correspond- 
ing condensation energy, if any, is many times lower than 
the condensation energy caused by the realistic electron- 
phonon interaction. 



II. ALTERNATIVE THEORY 

Based on the above analytical and numerical re- 
sults and numerous experimental facts supporting strong 
electron-lattice coupling [l5[ we trust that EPI should 
be the key in solving the high-T c problem. However, 
a quantitative analysis of the doping-dependent EPI has 
remained elusive in cuprates and related compounds. Re- 
cent observations of the quantum magnetic oscillations in 
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FIG. 2: Condensation energy per electron (in units of the 
hopping integral t) versus the amplitude of the superconduct- 
ing d-wave order-parameter in the 2D Hubbard model with 
the finite-range EPI for U/t — 8, and different EP coupling 
A. (Reproduced from Ref. [l3|, ©American Physical Society, 
2009.) 



trix element, M(q), and its magnitude are well known, 
\M(q)\ = j(q))tiw /V2N. Here 7(g) = y/ATTe 2 /Kflhuj a q 2 
is a dimensionless coupling, f2 is a unit cell volume, N 
is the number of unit cells in a crystal, luq is the opti- 
cal phonon frequency, and k = £ooeo/(eo — £00) (coo and 
eo are high-frequency and static dielectric constants in a 
parent polar insulator). 

Recently I have noticed that, in highly polarizable 
ionic lattices as cuprates with eo > 1 the bare long- 
range Coulomb and Frohlich interactions almost negate 
each other allowing for the analytical multi-polaron the- 
ory in the strong-coupling regime (2lj . The dielectric 
response function of strongly correlated electrons is apri- 
ori unknown. Hence one has to start here with a generic 
Hamiltonian including unscreened Coulomb and Frohlich 
interactions operating on the same scale since any ad-hoc 
assumption on their range and relative magnitude might 
fail, 



H = - ^(TijSsj + fi6ij)4cj + ^^2 



2 z —' £oolm - n| 

1^3 



hioofii [u(m, q)dq + H.c] + H : 



ph • 
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some underdoped [16| and overdoped [l7|] cuprate super- 
conductors has opened up the possibility for a quantita- 
tive assessment of EPI in these and related doped ionic 
lattices [llj]. The oscillations revealed cylindrical Fermi 
surfaces, enhanced effective masses of carriers (ranging 
from 2m e to 6m e ) and the astonishingly low Fermi en- 
ergy, which appears to be well below 40 meV in under- 
doped [3] and less than or about 400 meV in heavily 
overdoped cuprates jl7j . Such low Fermi energies make 
the Migdal-Eliashberg adiabatic approach to EPI inap- 
plicable in these compounds. Since carriers in cuprates 
are in the non-adiabatic (underdoped) or near-adiabatic 
(overdoped) regimes with their Fermi energy about the 
optical phonon energy, our multi polaron theory of su- 
perconductivity Q should be more appropriate here. 

The theory of dense polaronic systems in the interme- 
diate coupling regime has remained rather cumbersome, 
in particular, when EPI competes with strong electron 
correlations. Corresponding microscopic models with the 
on-site Hubbard repulsion and the short-range Holstein 
EPI have been studied using analytical and powerful nu- 
merical techniques [l5j . 

In most analytical and numerical studies both inter- 
actions have been introduced as input parameters not 
directly related to the material. Quantitative calcula- 
tions of the interaction matrix elements can be performed 
from pseudopotentials using the density functional the- 
ory (DFT) [13]. On the other hand, one can express 
the bare Coulomb repulsion and EPI through material 
parameters rather than computing them from first prin- 
ciples in many physically important cases 20] . In partic- 
ular, for a polar coupling to longitudinal optical phonons 
(the Frohlich EPI), which is the major EPI in polar 
crystals, both the momentum dependence of the ma- 



Here Tij = T(m — n) is the bare hopping integral, \x 
is the chemical potential, i ~ m, s and j = n, s' in- 
clude both site (m,n) and spin (s, s') states, tt(m, q) = 
(2A r )~ 1 / 2 7(g) exp(zq • m), Cj, d q are electron and phonon 
operators, respectively, hi — cjc, is a site occupation 
operator, and H p h = X) q ^o^q^q + 1/2) is the polar 
vibration energy. 

This Hamiltonian can be reduced to the polaronic "t- 
J p " Hamiltonian using two successive canonical transfor- 
mations [2l| . 

n = -^2(tijS ss > + fiSi^cjcj 

+ 2^2 J P ( m - n ) ( S m • S n + jh m h n ) , (4) 



'ms ' ss' i-ms' 



where S m = (1/2) £ 
erator (f are the Pauli matrices) 



is the spin 1/2 op- 

h m =J2 s ^h and A is 
the renormalised chemical potential. 

There is a striking difference between this polaronic t- 
J p Hamiltonian Eq.Q and the familiar t-J model derived 
from the repulsive Hubbard U Hamiltonian in the limit 
U ^> t omitting the so-called three-site hoppings and EPI 
[22I ]. The latter model acts in a projected Hilbert space 
constrained to no double occupancy. Within this stan- 
dard t-J model the bare transfer amplitude of electrons 
(t) sets the energy scale for incoherent transport, while 
the Heisenberg interaction (J oc t 2 /U) allows for spin 
flips leading to coherent hole motion with an effective 
bandwidth determined by J <C t. On the contrary in our 
polaronic t-3 p Hamiltonian, Eq.Q there is no constraint 
on the double on-site occupancy since the Coulomb re- 
pulsion is negated by the Frohlich EPI. The polaronic 
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hopping integral t leads to the coherent (bi)polaron band 
and the antifcrromagnetic exchange of purely phononic 
origin J p ^> t bounds polarons into small superlight inter- 
site bipolarons. Last but not least the difference is in 
the "+" sign in the last term of Eq.Q proportional to 
h m n n , which protects the ground superconducting state 
from the bipolaron clustering, in contrast with the "-" 
sign in the similar term of the standard t-J model, where 
the phase separation is expected at sufficiently large J. 

The reduction of Eq.(|3]) to Eq.(j4|) is based on the 
so-called "inverse coupling constant (1/A)" perturba- 
tion technique developed by us for the multi-polaron 
systems [3j, where the residual polaron-phonon interac- 
tion creates multi-phonon vertexes in the diagrammatic 
technique. Different from any model proposed so far, 
all quantities in the polaronic t-J p Hamiltonian (f?|) are 
defined through the material parameters, in particular 
tij = T(m — n) exp[— g 2 (m — n)] with 

2/ n 27re v— * 1 — cosfq ■ m) 

9 m = z ATO J2 T 1 5 

q 

and 



Here the high-frequency, and the static, eo dielectric 
constants as well as the optical phonon frequency, lvq and 
the bare hopping integrals in a rigid lattice, T(m) are 
measured and/or found using the first-principle Density 
Functional Theory [19( in a parent polar insulator. 

Solving the polaronic t- J p model with the appropriate 
material parameters, k,ujq and T(m) exhibits a phase 
transition to a superconducting state with the critical 
temperature in excess of 100K at optimum doping [2ll . 
[23| . The model provides un understanding of the spin 
and charge pseudogaps [23j and the unified parameter- 
free explanation of the observed oxygen-isotope efects on 
the critical temperature, the magnetic-field penetration 
depth, and on the normal-state pseudogap in underdoped 
cuprate superconductors [241 ] . 
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